Abstract. We propose a vacuum gap (VG) model which 
can be appHed uniformly for normal and high magnetic 
field pulsars. The model requires strong and non-dipolar 
surface magnetic field near the pulsar polar cap. We as- 
sume that the actual surface magnetic field Bg in pul- 
sars results from a superposition of global dipole field Bd 
and crust-anchored small scale magnetic anomaly Bm. We 
provide a numerical formalism for modelling such struc- 
tures of surface magnetic field and explore it within the 
' framework of VG model, which requires strong surface 
fields Bs> 10^^ G. Thus, in order to increase the resul- 
tant surface field to values exceeding 10^'^ G, in low mag- 
, netic field pulsars with Bd 10^'^ G it is required that 
^ ' Bm 3> Bd, with the same polarities (orientations) of B^ 
O , and Bm. However, if the polarities are opposite, the re- 
' sultant surface field can be lower than the dipolar surface 
. component inferred from the pulsar spin-down. Wc pro- 
P\J ' pose that high magnetic field pulsars (HBPs) with the in- 
ferred global dipole field Bd exceeding the so called photon 
splitting threshold Bcr ~ 4 x lO^'^ G, can generate observ- 
able radio emission 'against the odds', provided that the 
surface dipolar magnetic field Bd is reduced below Bcr by 
, the magnetic anomaly Bm of the right strength and po- 
T— I ' larity. We find that the effective reduction is possible if 
,; the values of Bd and Bm are of the same order of magni- 
' tude, which should be expected in HBPs with Bd > Bcr- 
The proposed VG model of radio emission from HBPs, in 
■ which pair production occurs right above the polar cap, is 
an alternative to the recently proposed lengthened space 



> 



Oh- 



O ■ charge limited fiow (SCLF) model, in which pair forma- 
^ , tion front is located at relatively high altitudes, where 
^ ' the dipole field is degraded below Be - Our model allows 
• • , high Bd radio-loud pulsars not only just above Bcr but 
^ ' even above 2 x 10"'^^ G, which is the upper limit for HBPs 
, within the lengthened SCLF model. 
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1. Introduction 



(1986| ) showed that in the crustal model it is only possible 



The properties of radio emission of typical pulsars strongly 
suggest that the magnetic field is purely dipolar, at least 
at altitudes r of several stellar radii R = 10^ cm, where the 



to form small scale surface field anomalies with a typical 



radio emiss ion is expected to originate (e.g. Kijak fc Gil 
1997 , 1998 , and references therein). However, this may 
not be a good description of the structure of the mag- 
netic field at the stellar surface. In fact, already Ruder- 
man & Sutherland (RS75) implicitly assumed that the ra- 
dius of curvature of field lines above the polar cap should 
be about 10^ cm, which is inconsistent with the global 
dipolar magnetic field. Several authors argued on theoret- 
ical grounds that the magnetic field could be produced 
by currents fiowing in thin crustal layers of the neutron 
star, which would generate non-dipolar fields at the sur 
face (e.g 
lik, 1991 



size of the order of 100 meters. Gil & Mitra (2001) demon- 
strated that such 'sunspot' like magnetic field structures 
on the polar cap surface help to sustain VG-driven radio 
emission of pulsars. Here we envisage the scenario where 
the magnetic field of neutron star is non-dipolar in na- 
ture as a superposition of the fossil field in the core and 
the crustal field structures. The crust gives rise to small 
scale anomalies which can be modelled by a number of 
crust anchored dipoles oriented in different directions (e.g. 



Blandford, Applegate & Hernquist, 1983; Arons, 1993) 



Blandford, Applegate & Hernquist, 1983; Kro 



Ruderman, 1991; Arons, 1993; Chen & Ruder- 



man, 1993; Geppert & Urpin, 1994; Mitra, Konar & Bhat 



tacharya, 1999). A lot of observational evidence has also 
been pfescn ted.' Pago & Sarmicmo (1990 ) and Bulik ci al 
(1992| , 1995 ) reported that imerpreration of their analysis 
of X-ray pulsars suggested small scale magnetic anoma- 
lies on the polar cap, which would strongly deviate the 
surface field from the purely dipolar configuration. It is 



The superposition of global dipole and local anomaly is 
illustrated in Fig. 1, where for clarity of presentation only 
one local, crust associated dipole is marked. 

Formation of dense electron-positron pair plasma is 
essential for pulsar radiation, especially (but not only) 
at radio wavelengths. Purely quantum process for mag- 
netic pair production 7 —> e~e ^ is commonly invoked 



as a source of this plasm a (e.g. Sturrock, 1971 ; Ruder 



man & Sutherland, 1975). However, at superstrong mag- 
netic fields close to the so-called quantum field Bq = 
4.4 X 10^'^ G, the process of free e~e+ pair production 
can be dominated by the phenomenon of photon split- 



b elieved that thermal X-rays from the polar cap surface 
are good diagnostic tool to infer the structure of the sur- 
face magnetic field. Several similar arguments in favour of 
non-dipolar nature of surface magnetic field can also be 



ting (Adler et al., 1970; Bialynicka & Bialynicki et al. 
1970| ; Baring & Harding, 1998) and/or bound positron- 
ium formation (Usov & Melrose, 1995, 1996| ). While the 



latter process can reduce the number of free pairs at mag- 



found in 


Becker & Triimper (1997|); 


Cheng, Gil & Zhang 


(1998); 1 


iudak & Dyks (1999); Cheng & Zhang (19991); 


Thompson & Duncan (1995|, |1996 


), and 
)• 


Murakami et al. 


(1999), and Tauris & Konar (2001 





netic fields B>0.lBq (e.g. [Baring fc Harding, 200l| ), the 



former one can entirely suppress the magnetic pair pro- 
duction at B> 10^^ G, provided that photons polarized 
both parallel and perpendicular to local magnetic field di- 



Wo ftjer (1964) proposed that the magnetic field is the 
fossil field of the progenitor star amplified during the col- 
lapse and anchored in the superfluid core of the neutron 
star. We will assume in this paper that the magnetic field 
was also generated in the outer crust during or shortly af- 
ter the neutron star was formed by some u nspecified mech- 
anism (e.g. thermomagnetic instabilities; Blandford, Ap- 



rection can split (e.g. Baring, 2001; Baring fc Harding 



2001). This assumption will be implicitly kept throughout 
this paper. Under these circumstances one can roughly de- 
fine a photon splitting critical line Bcr Bq and expect 
that there should be no radio pulsar above this line on 
the Bd- P diagram, where Bd = 6.4 x lO^^ (PP)!/^ q 
is the dipole surface magnetic field estimated at the pole 
from the pulsar p eriod P and its derivat ive P ( Shapiro fc 



plegatt fc Hernquist, 1983 ). Urpin, Levshakov fc lakovlev Tcukolsky, 1985 ; Usov fc Melrose, 1996 ). This death-line 
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is more illustrative than quantitative. In fact, a number 
of specific model dependent death-lines separating radio- 
loud from radio-quiet pulsars are available in the litera 
ture 



8 


5aring & Harding, 1998, 2001 


7 


Zhang & Harding, 


2000b) 




2001 


). All these slightly period dependent death- 



lines cluster around Bq on the Bd — P diagram, and hence 
the quantum field is conventionally treated as a threshold 
magnetic field above which pulsar radio emission ceases. 
In this paper we also use this terminology, bearing in mind 
that the photon splitting threshold realistically means a 
narrow range of magnetic fields around the critical quan- 
tum field -Bp ^ 4 X 10 ^^ G, certainly above 10^^ G (see 
review by [Baring, 2001] ). For convenience, in all numerical 
examples presented in Figs. 2-6 and subsequent discus- 
sions we assume the threshold magnetic field Bcr ~ Bq. 

In order to produce the necessary dense electron- 
positron plasma, high voltage accelerating region has to 
exist near the polar cap of pulsars. Two models of such ac- 
celeration regions were proposed: stationary space charge 
limi ted fiow (SCLF) models (^harleman, Arons fc Fawley, 



1978; Arons & Sharleman, 1979; Arons, 1981) in which 



charged particles fiow freely from the polar cap, and highly 
non-stationary vacuum gap (VG) models (Ruderman & 

198C| : 



and, 1975; Cheng & Ruderman, 1977 



Gil & 



Suther 

Mitra, 2001 ) in which the free outflow of charged particles 
from the polar cap surface is strongly impeded. In the VG 
models the charged particles accelerate within a height 
scale of about polar cap radius of ~ 10"* cm, due to high 
potential drop across the gap, while in the SCLF models 
particles accelerate within a height scale of a stellar ra- 
dius ~ 10^ cm, due to the potential drop resulting from 
the curvature of field lines and/or inertia of outstreaming 
particles. In both models the free e~e"'' pairs are created 
if the kinematic threshold e-f ■ sin 9t = 2mc^ is reached or 
exceeded and the local magnetic field is lower than the 
photon splitting threshold B ^ Bcr, where = huj is 
the photon energy and 9t is the propagation angle with 
respect to the direction of the local magnetic field. 

Recent discovery of high magnetic field pulsars (HBPs) 
however has challenged the existing pair creation theories. 
A few HBPs found the inferred surface dipolar fields above 
the photon splitting level: PSRs Jl 119-6127, J1814-1744 
and J1726— 3530 (Table 1). Moreover, yet another strong 



Table 1. Radio- loud HBPs with inferred magnetic field 
Bd = 6.4 X 10i^(PP)i/2 Q higher tha n critical quantum 
field Brr = 4.4 X 10^3 q (^fter Table 1 in [Zhang fc Harding 



source 




P(s) 


P (s/s) 




Bd (G) 


PSR J1814- 


1744 


3.98 


7.43 X 10" 


la 


1.1 X 10'* 


PSR J1119- 


6127 


0.41 


4.02 X 10" 


12 


8.2 X 10^^ 


PSR J1726- 


3530 


1.11 


1.22 X 10" 


12 


7.4 X 10^3 



erators (thus no e~e^ pair production), while the HBPs 
can, with a difference attributed to the relative orienta- 
tions of rotation and magnetic axes (neutron stars can be 
either parallel rotator (PRs) with HB^ > or antiparallel 
rotator (APRs) with JIB^ < 0, where Q is the pulsar spin 
axis and is the magnetic field at the pole). If the pho- 
ton splitting suppresses completely the pair production at 
the polar cap surface, then the VG inner accelerator can- 
not form, since the high potential drop cannot be screened 
at the top of the acceleration region. Hence, ZHOO argued 
that in high magnetic field regime (Bd > Bcr) the pair 
production process is possible only if the SCLF acceler- 
ator forms. In fact, such SCLF accelerators are typically 
quite long and their pair formation front (PFF) can oc- 
cur at high altitudes r, where the dipolar magnetic field 
Bd oc r^^ has degraded below the critical value Bcr- Fur- 
thermore, ZHOO demonstrated that such lengthened SCLF 
accelerator in magnetar environment can form only for 
PRs and not for APRs. Consequently they concluded that 
the radio-loud HBPs are PRs with developed lengthened 
SCLF accelerator, while the radio-quiet magnetars (AXPs 
and SGRs) represent APRs. It is worth emphasizing here 
that ZHOO developed their model under the assumption 
that the magnetic field at the surface of HBPs is purely 
dipolar. 

In this paper we propose an alternative model for 
radio-loud HBPs based on highly non-dipolar surface mag- 
netic field, in which the photon splitting within the VG in- 
ner acceleration region does not operate even if the dipole 
magnetic field exceeds the critical value Bcr at the po- 
lar cap. Thus our model requires that HBPs are APRs, 



field nc futron star PSR J1846-0258 with Pd - 5 x 10^^ G 
was discovered (potthelf et al., 200C| ), which seems to 



which is a consequence of the VG scenario (e.g. Ruder 



be radio-quiet ( Kaspi et al., 1996[ ), although its X-ray 



man fc Sutherland, 1975; Gil fc Mitra, 2001). This model 



emission is apparently driven by dense e e"*" pair plasma 



(e.g. Cordes, 2001). However, one should keep in mind 
that the actual threshold due to photon splitting and/or 
bound positronium formation can be well below the criti- 
cal field Bcr ^ 4x10^"^ G, indicating that all high magnetic 
field radio pulsars with Bd > 10^^ G pose a challenge. 
To evade the photon splitting problem for these pulsars 
Zhang fc Harding (2000a, ZHOO hereafter) proposed "a 
unified picture for HBPs and magnetars" . They argued 
that radio-quiet magnetars cannot have active inner accel- 



is a follow-up work of Gil & Mitra (2001), who argued 
that the VG can form if the actual surface magnetic field 
is about 10"^'^ Gauss. In other words, they assumed that all 
VG-driven radio pulsars have a very strong, highly non- 
dipolar surface magnetic field, with a strength more or less 
independent of the value of the global dipole field inferred 
from the magnetic breaking law. Thus, if Bd <^ 10^^ G 
then Bs > Bd and if Bd> 10" G then B^ Bd\ how- 
ever, in any case Bg < Bcr ~ Bq. We argue that such 
strong surface field anomalies can increase low dipolar field 
Bd ^ Bcr in normal pulsars to values exceeding 10^^ G 
(required by conditions for VG formation - see Gil fc Mi- 
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Fig. 1. Superposition of the star centered global magnetic 
dipole d and crust anchored local dipole m placed at = 
{vs ^ -R, = 9r) and inclined to the z-axis by an angle 
9m- The actual surface magnetic field at radius vector r — 
{r, 9) is B, = Bd+Bm, where = 2d/r^, B,„ = 2m/|r- 
Tsl, r is the radius (altitude) and 9 is the polar angle 
(magnetic colatitude) . R is the radius of the neutron star 
and L is the crust thickness. 

tra 2001) if global and local surface fields have the same 
polarities, or reduce very high dipolar field Bd> Bcr in 
HBPs if both these components are of comparable values 
and have opposite polarities. 

2. Modelling the Surface Magnetic Field 

We model the actual surface magnetic field by superposi- 
tion of the star centered global dipole d and a crust an- 
chored dipole moment m, whose influence results in small 
scale deviations of surface magnetic field from the global 
dipole. According to general situation presented in Fig. ^ 

Bs = Bd -I- Bm, (1) 

where 

B.= f^, ^,oV (2) 



r and 9 are star-centered polar co-ordinates, and 



B„ 



3 (r - r^) (m- (r - r^)) - m |r - 



(3) 



The global magnetic dipole moment d = {l/2)BdB? , 
where Ed = 6.4 • 10^^ (P • P)^/^ G is the dipole component 
at the pole derived from the pulsar spin-down rate, and the 
crust anchored local dipole moment m = (l/2)i?„iAi?^, 
where Ai? ^ 0.05i? is the characteristic crust dimension 
[R = 10^ cm). We use spherical coordinates with z axis di- 
rected along the global magnetic dipole moment d. Thus, 

= {rs,9r, (j)r) and m = (m, 6'™, </)„). 

To obtain the equation of the open magnetic field lines 
we first define the boundary of the open field lines at an 
altitude where the magnetic field should be a pure dipole 
and we have chosen as a starting altitude r = 5QR. Then 
we solve the system of differential equations 
d9 _ Bl + B^ 

^-Z71^m^ = 0u (4) 
dr 



r {B^ + Bl: 
BT 



(5) 



r [B'j: + P™) sin 6* 
with the initial conditions B,,,, = defined &t r — 50R 
(already at r — 5R the ratio B„i/Bd ^ lO"**) and trace 
the field lines down to the stellar surface {r — R). Here 

P™ = 



1 



Be - - 

and 
D 



£,2.5 



D 



1 



D2.5 



(3Tr^ - 3Tr + Drrir) 
Dme) , 



(6) 



cost 



— 2rsr (sin 9r sin 9 cos {(j) — (pr) + cos t 
T = rrirr - (m^^'r + i^erl + m^r^) , 

~ Ts (sin6'r sinfl cos (0 — 1/)^) -f cos^^ cos 9) 
rg = Ts (sin 9r cos 9 cos {(j) — (pr) — cos 9r sin 9) 

= — sin^r sin (0 — (/)r) 
rrir — m {sin 9m sin 9 cos (0 — (pm) + cos 9m cos I 
mg — m (sin6',„ cos 9 cos {(f> — (pm) — cos 9m sin( 
= -TO sin 6*™ sin {4> - 4>m) ■ 
The curvature pc = of the field lines (where 5R 
is the radius of curvature presented for various cases in 
Fig.H) is calculated as 



(7) 



ds 



d^T ds dvd^s 
dr"^ dr dr dr"^ 



p,^{S,)-^Jl + Jl + Jlf\ 
where 

J = ^(f + r2e2+r2$2 sin2 Q) , 

Jl = — X1S2, J2 = ^2'5'l ~ Y1S2, 

J3 = Z2S1 — ^15*2 

Xi = sin cos (/> + r0i COS 6* COS (/) — r$i sin0 sin (/), 

Yi — sin(?sin(/) + r0i cos6'sin0 + r"I>i sin^cosc/), 

Zi = cos9 — rQisin9, 

X2 = (281 -f're2)cos6'coS(?!)- (2$i -f r$2)sin6'sin(: 
- r (61 + $1) sin cost/) - 2r6i$i cos 6* sin 0, 



(8) 
(9) 
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Y2 = 

Z2 = 

51 = 

52 = 



(29i + r92) cos0sin0 + (2<i>i + r$2) sin0 cos < 
r (Oj + $1) sin 6* sin + 2rQi^i cos 6* cos 0, 
—01 sin 9 — Qi sin 9 — r02 sin 9 — rOl cos 9, 



m=-10 d 



1 

Si' 



r'^^l sin 9, 
26182 +r$2 sin^o 



92 = 



+ r2$i$2 sin^ 6* + r^Oi^i sin6lcos6') . 
dOi d$i , , 

For simplicity, in this paper we mostly consider an 
axially symmetric case in which both d and m are di- 
rected along z-axis (parallel or antiparallel), thus 9^ — 
0m = (f'r = 4'm. = 0. Also, for Convenience, m is ex- 
pressed in units of d. We use normalized units in which 
d = P = R — 1 and — 0.95. (see caption of Fig. 2 for the 
normalization convention). All calculations are carried out 
in three-dimensions, although, for clarity of graphic pre- 
sentation, in figures we present only two-dimensional cuts 
of the open field line regions. 

2.1. High magnetic field pulsars 

As mentioned above, the formation of VG inner acceler- 
ator requires very high magnetic field B> 10^^ G on the 

T996HGil 



surface of the polar cap (Usov & Melrose, 1995 



fe Mitra, 2001). This can be achieved not only in pulsars 
with high dipolar field Bj^> 10^^ G. In fact, some of the 
low field pulsars with B^i ^ 10^^ G can have surface field 
Bs> 10^3 G if B„ > Bd (thus m > 1.25 x lO^* d). We 
discuss such normal, low field pulsars later in this paper. 
Presently let us consider the HBP with a dipolar surface 
field at the pole Bd = 6.4 x 10^^{P ■ P)i/2 G exceeding 
the photon splitting limit Bcr ~ Bg. If all photon split- 
ting modes operate, such pulsar should be radio-quiet. Al- 
ternatively these pulsars could be radio-loud if the effec- 
tive surface field is reduced below Bcr- Such scenario can 
be achieved if the polarities of magnetic moments d and 
m are opposite, that is d and m are antiparallel. Fig. || 
presents a case with m = — 10~^d and AR/R = 0.05. The 
actual surface values of B"^ = y/ {B^Y + {B^Y as well as 
radial components of B;' = ■ R/i? and i?^ = ■ R/i? 
are presented in the lower panel of Fig. || (note that all 
radial components are positive and that the total i?** is al- 
most equal B^ in this case) . At the pole (radius r = R and 
polar angle 6* = 0) the ratio B^/Ba = (m/d) ■ (R/ARf = 
0.8 and thus Bs =Bd-B^ = 5^(1-0.8) = 0.25^- As one 
can see from this figure, all surface field lines between —9s 
and +9s are open, but the ratio Bg/Bd increases towards 
the polar cap edge, reaching the value of about 0.5 in the 
region between polar angles \9d\ and \9s\. The ratio Bm/Bd 
is also about 0.5 in this region. Thus, the global dipolar 
field {Bd = 2 in our units) is effectively reduced between 2 
and 5 times in different parts of the polar cap (defined as 
the surface area from which the open magnetic field lines 
emanate). This means that the ratio B^^/Bd ranges from 




Fig. 2. Structure of surface magnetic field for a super- 
position of the global star centered dipolar moment d 
and crust anchored dipole moment m = — 10~^d (Fig. |l|). 
The open dipolar field lines (solid) and the actual surface 
open field lines (dashed) are shown in the upper panel. 
The horizontal axis is labelled by an azimuthal angle 9 
(magnetic colatitude), which measures the polar cap ra- 
dius. For purely dipolar field lines the polar cap radius 
rd ~ R ■ sm9d, which for pulsar period P = 1 s is about 
0.014 radians (thus rd ~ 1.4 • 10^ cm). The actual polar 
cap is broader with the last open lines emanating at the 
polar angles 9s « 0.023 (thus the actual polar cap radius 
Vs — 2.3 • 10^ cm — 1.657"^). The actual open surface field 
lines (solid) reconnect with dipolar ones (dashed) at dis- 
tances z — (r/R) ■ cos 9 w 1.2, where r is the radius and 
9 < 0.025 radians. In the lower panel the surface values 
(r = R) of both dipolar field (dashed horizontal) and the 
actual field (solid line) are shown. The radial components 
B^ = 2dcos9/R^ 2 and B^ = {Bd + B„) • r/r, and 
total values Bs = ^y{BfY + {BfY are presented (where 
c? = i?, = 1 is assumed for convenience) . 
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0.5 to 0.8 across the polar cap. The actual polar cap is 
broader than the canonical dipolar polar cap (two dashed 
vertical lines correspond to last open dipolar field lines 
emanating at the polar angles 9d = ±0.014 radians for 
typical period P — Is). The ratio of actual to dipolar po- 
lar cap radii is 6s /6d ^ 5/3 in this case. Thus, using the 
argument of magnetic flux conservation of the open field 
lines, one can say that an effective surface magnetic field 
of the polar cap is about 2.8 times lower than the dipolar 
surface field measured from the values of P and P. If the 
estimated dipole field « 10^^ G (like in the case of 
PSR J1814-1744, Table 1) then the actual surface field 
at the pole is only Bg ^ 2.5 • 10^'^ G, well below the pho- 
ton splitting death line Be- — 4.4 x 10^^ G. Such pulsar 
can be radio-loud without invoking the lengt hened SCLF 
accelerator proposed by ZHOO. As shown by Gil & Mitra 
(200l|)] in such strong surface magnetic field the vacuum 
gap accelerator can form, which implies low altitude co- 
herent radio emission ( [Melikidzc, Gil fc Pataraya, 2000| , 



m=-2x10 d 



see section 3 in this paper) at altitudes rem ~ 50i? (for 
a typical pulsar with P 1 s) in agreement with obser- 



vational constraints on radio emission altitudes (Cordes 
1978|, |]|992|; |Kijak fc Gil, 1997|, |1998|; |Kijak, 20011). 



Fig. ^ presents another case of opposite polarities 
m = — 2 X 10~^d, with a magnitude of m two times 
stronger than in the previous case (Fig. ||). Again for 
AR/R - 0.05, B,n/Bd = {m/d){R/ARf = 1.6 at the 
pole (r = i? and 61 0) and B^ = Ba - B,n = -O.GBd- 
The negative sign of the ratio Bg/Bd means that the sur- 
face magnetic field Bs is directed opposite to near the 
pole, that is the circumpolar field lines with polar angles 
-9si < 9 < Osi (where 9si ~ 0.031) are closed. The last 
open surface field lines (solid) emanating at polar angles 
9s = ±^s2 (where 6*32 ^ 0.037) reconnect with last open 
dipolar field lines (dashed) at altitudes 2; ~ 1.2 (thus about 
2 km above the surface). The actual polar cap, which is the 
surface through which the open magnetic field lines em- 
anate, has a shape of a ring (0.031< |^5|< 0.037) located 
outside the circle of dipolar polar cap with angular ra- 
dius 9d = 0.014 (or diameter ft; 6*^ • i? « 1.4 • lO^* cm). 
Again, the magnetic flux conservation argument leads to 
Bs/Bd = {9X-9\)/9l = (0.0372-0.03l2)/0.0142 = 0.48, 
thus Bs is about 0.55^ within the ring of the open 
field lines (thus \Bm\IBd ^ 1.6 in this region. The ac- 
tual values of surface magnetic field (radial B^ and to- 
tal = {BrY + i^Y) ^'"S shown as solid lines in 
the lower panel of Fig. 0, in comparison with radial com- 
ponents of dipolar field Bd (dashed horizontal line). As 
one can see, B^ < Bf = 2m cos 9 /R^ and B^ is neg- 
ative for \9\ < 0.022. If the dipolar surface component 
of pulsar magnetic field ~ 7 x 10^^ G (hke in PSR 
J1726— 3530, Table 1), then the actual surface magnetic 
field B 
old. 




4 X 10^^ G, below Lhc phuLuii spliLLliig Lhiusli- — Mitra 2001) 



Fig. 3. As in Fig. | but for m = -2 x lO^^'d. 



2.2. Normal pulsars 

Fig. ^ presents a case with m = 2 x 10~'*d in which both 
magnetic moments have the same polarity. Obviously in 
such case, the surface magnetic field will be stronger as 
compared with pure dipole (m = 0). Flux conservation ar- 
gument gives surface magnetic field Bs/Bd — {9d/9sY = 
(0.014/0.008)2 - 3 (the ration B.,n/Bd< 2 x 10^^8000 = 
1.6). Thus the actual surface field is about 3 times stronger 
than the inferred dipolar field Bp = 6.4 x 10^^{P-PY/^ G. 
It is interesting to compare this case with the previous one 
(m = — 2 X 10~^d presented in Fig. |^), in which the actual 
surface field Bs is about 2 times weaker than the global 
dipolar surface field at the polar cap. Such cases of in- 
creasing an effective magnetic field can be important in 
normal pulsars with low dipolar field Bd ^ 10^'^ G (Gil & 



It is then interesting to examine how different polar- 
ities of d and m would influence normal pulsars with 
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m=2x10 d 



m=-4x10 d 




1.14 




Fig. 4. As in Fig. | but for m = 2 x lO^^d. 



Fig. 5. As in Fig. | but for m = -4 x lO^^'d. 



Bd < B„. If m/d ~ (AR/R)^ thus S,„ - Bd then of 
course Bg can be slightly lower then Bd, as in the case 
of HBPs (Fig. In such case, however, the VG can- 
not form. In fact, as arg ued by |Gil fc Mitra (200l| ), the 
formation of VG requires that Bg is close to 10^'^ G or 
even above, thus Bm ^ Bd is required in normal pul- 
sars (see also Gil et al. 2001). Fig. ^ illustrates a case 
of high surface magnetic field with Bm ^ Bd, in which 
VG can apparently form. As one can see from this fig- 
ure, the values of Bs at the ring-shaped polar cap are 
close to dipolar values Bg = Bd{r = R, \9\ ^ 0.05). 
One can show that this is a general situation, that is 
Bs ~ Bd no matter how much Bm exceeds Bd at the 
pole. This follows from the fact that the angular location 
9 of the polar cap ring increases with the increasing ratio 
Bm/Bd w {m/d){AR/Rf ~ 8 x 10'^[m/d). For example, 
in the case presented in Fig. |^ m = —4 x 10^'^d (thus for 
Ai? = 0.05i? we have Bm ~ iBd at the pole) and the last 
open field lines emanate at polar angles Og ~ ±0.055 radi- 



ans, or at polar cap radii i?p ~ 6 x 10** cm (for _P = 1 s). 
Thus, the narrow polar cap ring is located far from the 
local dipole m, whose influence is weak at this distance. 
The circumpolar field lines between polar angles —0.053 
to -1-0.053 are closed. 

Thus, we conclude that the actual pulsar surface mag- 
netic field Bg can significantly differ (say by an order of 
magnitude) from the inferred dipolar field Bd only in the 
case when the polarities of the global d and local m dipole 
(Fig. ^ are the same, as illustrated in Fig. |4| If this is the 
case, then Bg can largely exceed Bd, which seems to be 
important from the viewpoint of vacuum gap formation 
requiring B> G (see pil fc Mitra, 200l| ). Therefore, 
in normal VG driven radio pulsars the polar cap should be 
circular, or at least filled - if the axial symmetry does not 
hold. The ring-shaped polar cap can occur only in nor- 
mal pulsars with Bd< B^r and in radio-loud HBPs with 

Bd> Bcr- 
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Gil et al. (2001) explored consequences of the vac- 
uum gap model interpretation for drifting subpulses ob- 
served in PSR B0943-I-10, in which 20 sparks move cir- 
cumferentially around the perimeter of the polar cap, each 



completing one circulation in 37 pulsar periods (Desh- 
pande |fc Rankin, 1999| , |2GGl| ). pil et al. (200lD consid- 



ered both the curvature radiation (CR) and resonant in- 
verse Compton radiation (ICS), seed photons as sources 
of electron-positron pairs and determined the parame- 
ter space for the surface magnetic field structure in each 
case. For the CR-VG the surface magnetic field strength 
Bs > 2 X 10^'^ G and the radius of curvature of surface field 
lines 0.6 xlO^cm < TZ < 1.2xl0^cm, while for the resonant 
ICS-VG Bs>2x 1013 G and lO^cm < 7^ < 3 x lO^cm (of 
course, in both cased Bs < Bg ^ 4.4 x lO^^ q) xhc CR- 
VG with such curved surface magnetic field does not seem 
likely (although it cannot be excluded), while the ICS-VG 
gap supported by the magnetic field structure determined 
by the parameter space determined above guarantees a 



uum gap formation and drifting subpulse phenomenon. 
We model the actual surface magnetic field as the super- 
position of the global star-centered dipole and local crust- 
anchored dipoles Bg = -I- J2i Bmi ~ B^ -I- B„io, where 
B„io is the local dipole nearest to the polar cap centre 
(Fig. 1). Such model is quite general, as it describes the 
magnetic field structure even if the star-centered dipole is 
negliglible at the star surface. In such a case the surface 
dipole field 'B^ (inferred from P and P measurements) is 
a superposition of all crust-anchored dipoles calculated at 
far distance and projected down to the polar cap surface 
according to the dipolar law. 

We propose a model for radio-loud HBPs with high in- 
ferred dipolar magnetic field Bd > 10^^ G, even exceeding 
the critical value Bcr ^ 4 x lO^^^ G. Given the difficulty 
that at strong magnetic field the magnetic pair creation 
process is largely suppressed, the puzzling issue remains 
how these HBPs produce their e~e~^ pair plasma necessary 



system of 20 sparks circulating around the perimeter of 



for generation of the observable radio emission. Zhang & 
Harding (2000a) proposed a "lengthened version" of the 



the polar cap by means of the E x B drift in about 37 stationary SCLF model of inner accelerator (e.g. Arons fc 



pulsar periods. 



Further Gil et al. (2001 ) modelled the magnetic field 
structure determined by the ICS-VG parameter space 
(specified above), using the numerical formalism devel- 
oped in this paper. Since Bd = 6.4 x 10^^{P ■ P^/'^G = 
4 X 10^2 Q in this case, then to obtain B^ - (2-=-3) x 10" G 
one needs Bm ^ Bd and the same polarity of both com- 
ponents. Following the symmetry suggested by the ob- 
served patterns of drifting subpulses in PSR B0943-I-10, 
the local dipole axis was placed at the polar cap center. 
A number of model solutions corresponding to rg 0.97 
and m - (1 ^ 2) x 10"^ d and satisfying the ICS-VG 
parameter space, was then obtained. As a result of this 
specific modelling pil et al. (2001 ) obtained a number of 
interesting and important conclusions: (i) The conditions 
for the formation of the ICS-VG are satisfied only at pe- 
ripheral ring-like region of the polar cap, which can just 
accommodate a system of 20 E x B drifting sparks, (ii) 
The surface magnetic field lines within the actual gap are 
converging, which stabilizes the E x B drifting sparks by 
preventing them from rushing towards the pole (as op- 



Sharleman, 1979), in which the pair formation front De- 



posed to the case of diver ging dipolar field (e.g. Fillipcnko 
& Rac hakrishnan, 1982 ). (iii) No model solutions with 
~ (3 4) X 10^3 G and 7^ - (0.6 ~ 1.2) x 10^ cm, 
could be obtained which corresponding to the CR-VG pa- 
rameter space, which in turn favors the ICS-VG in PSR 
B0943-I-10. 



3. Discussion and conclusions 

We argue in this paper that a putative presence of strong 
non-dipolar magnetic field on the neutron star surface can 
help to understand recently discovered radio pulsars with 
magnetic field above the photon splitting threshold, as 
well as to understand long standing problems of the vac- 



curs at altitudes r high enough above the polar cap that 
Bd ^ Bcr{R/r)^ degrades below Bcr, thus evading the 
photon splitting threshold. Our VG model is an alterna- 
tive to the lengthened SCLF model, with pair creation 
occurring right at the polar cap surface, even if magnetic 
field exceeds Bcr- We have assumed that the open surface 
magnetic field lines result in an actual pulsar from super- 
position of the star centered global dipole moment and a 
crust anchored local dipole moment. We argued that if the 
polarities of these two components are opposite, and their 
values are comparable, then the actual value of the surface 
magnetic field Bs can be lower than the critical field Bcr , 
even if the global dipole field Bd exceeds the critical value. 
Thus, the creation of electron-positron plasma is possi- 
ble at least over a part of the polar cap and these high 
magnetic field neutron stars can be radio- loud (HBPs). In 
fact, one should expect that in HBPs, in which by defi- 
nition Bd> Bcr ~ = 4.4 X 10" G, the ratio Bm/Bd 
should be of the order of unity, since Bg = Bd + Bm and 
10" G< Bs< 4 X 10" G. 

Within our simple model of non-dipolar surface mag- 
netic field Bs one should expect that both cases m ■ d > 
and m • d < will occur with approximately equal prob- 
ability. However, from the viewpoint of observable radio 
emission only the latter case is interesting in HBPs with 
Bd> Bcr- In fact, when m-d > then the surface magnetic 
field Bs > Bd ^ Bcr (Fig. ^ and the photon splitting 
level is highly exceeded. For m • d < we have two pos- 
sibilities: (i) if m/d< (AR/R)^ thus Bm< Bd at the pole 
{r = R,6 = 0) then the polar cap (locus of the open 
field fines) is circular (Fig. |); (ii) if m/d > {AR/Rf 
thus B„i > Bd then part of the circumpolar field lines 
are closed and the actual polar cap has the shape of ring 
(Fig. H). In both above cases (i) and (ii), the actual sur- 
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face magnetic field Bs at the polar cap (or at least part 
of it) can be lower than Bcr, even if Bd exceeds Bcr- The 
values of Bm and Bd> Bcr should be comparable to make 
reduction of strong surface field Bg below B^ possible. 
In our illustrative examples presented in Figs. 2 and 3 
(corresponding to the same pulsar with P = 1 s and 
Bd = 2dR-^ = 6.4 X 10^^ (F • F)^/^ ^^^^ ^.^tios 

Bm/Bd ranging from 0.5 to 1.6. These values could be 
slightly different, say by a factor of few, thus we can say 
that the ratio Bm/Bd should be of the order of unity. If 
Bm/Bd » 1, then the reduction of surface dipole field 
is not effective (see example presented and discussed in 
Fig. 5). On the other hand, the case with Bm <Si Bd 
is not interesting, as it represents a weak surface mag- 
netic field anomaly. Thus, among a putative population 
of neutron stars with Bd> Bcr, only those with the ratio 
Bm/Bd = {m/d){R/ AR)^ of the order of unity, and with 
magnetic moment m and d (Fig. 1) antiparallel at the 
polar cap surface, that is m • d < 0, can be detected as 
HBPs. Other neutron stars from this population of high 
magnetic dipole field objects should be radio-quiet. This 
probably explains why there are so few HBPs detected. 

Within the lengthened SCLF model there is an up- 
per limit around Bd = 2 x 10^"' G for radio-loud HBPs 
(ZHOO, Zhang 2001) . As ZHOO argued, detecting a pulsar 
above this limit would strongly imply that only one mode 
of photon splitting occurs. Without the alternative model 
of HBPs proposed in this paper, such detection would re- 
ally have great importance for the fundamental physics 
of the photon splitting phenomenon. In our VG based 
model there is no natural upper limit for the radio-loud 
HBPs. However, it is known that due to the magnetic 
pressure the neutron star surface would tend to 'crack', 
which should occur at magnetic field strengths approach- 
ing 10^^ G ( Thompson fc Duncan, 199?: ). It is unclear how 
the radio emission would be affected due to such cracking 
process. 

To illustrate the above argument, let us consider Fig. || 
which presents yet another case of opposite polarities 
m = -1.25 X lO-'id. With AR/R - 0.05 this gives 
Bm/Bd ~ 1-0 and Bg = Bd — Bm = at the pole 
(r = i?, 6* = 0). The dashed horizontal line at B = 0.2 
in the lower panel corresponds to the surface magnetic 
field Bs which is 10 times weaker than the global dipole 
compo nent Bd = 2 (not shown in the figure). Thus if, 



-1.25x10 d 




Fig. 6. As in Fig. | but for m = -1.25 x 10"*d. See also 
text for explanation. 



cumpolar area with \9\ < 6y the surface field region Bs can 
be even more than 10 times weaker than Bd- Thus, within 
our model one can expect a radio-loud HBP with Bd even 
exceeding 4 x 10^'* G. However, their radio-beams should 
be much narrower than those expected in normal pulsars, 
at least few to several times less than {rem/ RY^J^P~2^ 
degrees (w here Vcm is the radio emission altitude; Kijak & 



for ex ample, Bd — 4: x 10^^ G (well above the length- Gil, 1997| , |1998| ). This would make such sources difScult 



ened SCLF limit Bd — 2 x 10^^ G; such pulsar was not 
observed so far), then the actual surface field Bs is well 
below Bcr ~ 4 X 10^'^ G, at least at the inner part of the 
polar cap between ±9^ = 0.0035 rad. This "pair form- 
ing effective" polar cap is about 2.5 times smaller than 



the canonical polar cap with radius 6d = 0.014 rad, and tra, 2001 



to detect 

The dotted horizontal line at i? = 0.05 in Fig. M cor- 
responds to Bs = 10^3 Q foj. adopted = 4 x W'^ G. 
This value of the surface magnetic field is believed to be 
about the lower limit for VG formation (see pil 



Mi- 



about 7 times smaller than the entire polar cap with radius 
9s — 0.027 rad. Near the last open field lines at polar an- 
gles 0.027> \9\> 0.014 the actual surface magnetic field Bs 
is only about 2 times lower than Bd, while in a narrow cir- 



Gil et al., 2001). Thus, the shadowed area in 
Fig. 1^ represent a narrow hollow-cone above which the VG 
driven radio emission cannot occur. Similar hollow-cone is 
marked in Fig. which presents the case similar to that 
illustrated in Fig. ^, except the local dipole is shifted off 
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-1.25x10'''d; Ae = 0.005 




-0.02 -0.01 

-6. 



0.01 0.02 0.03 



Fig. 7, 

by A0 



As in Fig. ^ but with local dipolc shifted off center 
= 0.005 radians. See also text for explanation. 



center by A9 — 0.005 radians (corresponding to about 0.2 
of the actual polar cap radius). The dashed horizontal line 
at B — 0.4 corresponds to Bs = 4 x 10^^ G and the dotted 
horizontal line at i? = 0.1 corresponds to Bg = 10^^ G, 
both calculated for adopted Bd — 2 x 10^^ G. The polar 
angles 6vi and 9y2 correspond to —9v and +6v in Fig. ^, re- 
spectively. The Fig. ^ demonstrates that conclusions of our 
paper do not depend on where the local dipole is placed. 

The above arguments strengthen the possibility that 
some magnetars can also emit observable radio emission 
( ICamiloetal., 200C| ; ^hang fc Harding, 2000b"| ). It is there- 



fore interesting to comment on the apparent proximity of 
HBP PSRJ 1814-1744 (with Bp = 1.1 x lO^"* G) and AXP 
IE 2259-f586 (with Bp = 1.2 x lO^"* G) on the P - P di- 
agram. In both these cases the inferred surface magnetic 
field well exceeds the critical value Per ■ Within our model, 
the former object can be radio-loud if the strong local 
dipole has the opposite polarity with respect to the global 



-1 — ' — I — ' — [- 




1. m = 

2. m = -lo'd 

3. m = -2x10'd 

4. m = 2x10"'d 

5. m = -4x10"'d 

6. m = -1.25x10*d 



Fig. 8. Radius of curvature pc (in units of i? = 10^ cm) 
as a function of normalized altitude z = (r/R) ■ cos 9 for 
actual surface magnetic field lines corresponding to four 
cases presented in Fig. 1 1, |, | and |. For comparison, the 
radius of curvature of purely dipolar field lines (in pulsar 
with P=l s) is shown (line 1). 



one. The radio quiescence of the latter object can be nat- 
urally explained if the local dipole is not able to decrease 
the inferred dipole magnetic field below the photon split- 
ting death-line. Thus, either the polarities are the same or 
they are opposite but the local dipole is not strong enough 
to reduce the dipole surface field below Bcr- 

In Fig. H we show the radii of curvature of actual 
surface field lines compared with those of purely dipo- 
lar field (line 1) as a function of normalized altitude 
z — {r/R)cos9 above the polar cap. Within the polar 
gap at z < 1.01 (within about 100 meters from the sur- 
face) the curvature radii for all cases presented in Fig. |[ 
^, ^, ^ an d |6| have values of the order of few hun dred me- 
ters (see Urpin, Levshakov & lakovlev , 1986), suitable 



for curvature radiation driven magnetic pair production 
{pc — 1/5R < 10® cm, where K is the curvature of field 
lines). 

All model calculations performed in this paper corre- 
spond to the axisymmetric case in which one local dipole 
is placed at the polar cap center (except the case presented 
in Fig. 0) . In the forthcoming paper we will consider a gen- 
eral, non-axisymmetric case, including more local dipoles, 
each with different orientation with respect to the global 
dipole. Although this generalization will give more realis- 
tic picture of an actual surface magnetic field, it will not 
change our conclusions obtained in this paper. 

It should be finally emphasized that although the 
lengthened SCLF model for HHPs (ZHOO) can solve the 
problem of pair creation in pulsars with surface dipole 
field exceeding the photon splitting threshold, it does not 
automatically warrants generation of the coherent radio 
emission of such HBPs. The problem is that unlike in 
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the non-stationary VG model, where the low altitude ra- 
dio emission can be generated by means of two-stream 



instabilities (Asseo fc Mchkidze, 199g; Melikidzc, Gil fc 



Patara^a, 2000), the stationary SCLF inner accelerator is 
associated with the high altitude relativistic maser radi- 



ation (e.g. Kazbegi, MachabeU fc Mehkidze, 1991, 1992; 



Kazbegi et al., 1996). This radiation requires relatively 
low Lorentz factors 7p ^ 5 ~ 10 of a dense secondary 
plasma (e.g. Machabeli & Usov, 1989). It is not clear 



if such plasma can be produced within the lengthened 
SCLF accelerator with delayed pair formation taking place 
in a purely dipolar magnetic field, either by cu rvature 



radiation or by inverse Compton scattering (e.g. Zhang 
& Har lling, 200"0b| ) processes. Moreover, the relativistic 



maser coherent radio emission requires a relatively weak 
magnetic field in the generation region. With the surface 
dipole field Bd ~ 10^'' G, such low field may not exist at 
reasonable altitudes (about 50% of the light cylinder ra- 
dius Rl = cP/2tt) required by the physics of correspond- 



ing instabilities (|Kazbcgi 


Machabeh & Melikidze, 1991 


1992; 


Kazbegi et al., 1996 


). Thus, if one assimrcs that the 



radio-loud HBPs are driven by the SCLF lengthened ac- 
celerator as proposed by ZHOO, they might not be able to 
generate observable coherent radio emission. This contra- 
diction seems to be a challenge for the lengthened SCLF 
scenario for HBPs. In our VG based model the low alti- 
tude (r ^ Rl) radio emission of HBPs is driven by just 
the same mechanism as the one most probably operating 
in typical radio pulsars (e.g. soliton curvature radiation 
proposed recently by Melikidze, Gil fc Pataraya, 2000| ) . In 
fact, the HBPs show apparently normal radio emission, 
with all properties typica l for characteristic pulsar radia- 
tion (ICamilo et al., 2000|) . 
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